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It is well known that supersymmetry (SUSY) gives neutron and electron electric 
dipole moments {dn and de) which are too large by about 10^. If we assume a SUSY 
model cannot contain fine-tunings or large mass scales, then one must require that 
the SUSY breaking mechanism give real soft breaking parameters, in which case the 
minimal SUSY model has no CP violation other than from the CKM matrix (besides 
possible strong CP violating effects). We show that in non-minimal SUSY models, 
a moderate amount of CP violation can be induced through one loop corrections 
to the scalar potential, giving an effective phase of order 10~^, and thus implying 
dn and d^ can be near their current experimental bounds naturally. This moderate 
amount of SUSY CP violation could also prove important for models of electroweak 
baryogenesis. We illustrate our results with a specific model. 
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1 Introduction 



Predictions for CP violating effects in supersymmetric (SUSY) theories have often 
been discussed with a certain ambiguity. On the one hand, it is well known that 
when the complex quantities in the theory are allowed to have phases of order unity, 
the predicted neutron and electron electric dipole moments {dn and de) are typically 
too large by perhaps 10^ [|T], ^ ^ |, |^. In order to avoid this, the relevant quantities 
are often chosen to be real, in which case the theory predicts no non-Standard Model 
CP violation (CPV) and negligible dn and de f^. On the other hand, it has often 
been assumed that the observation of dn around the current limit of 10~^^ecm 
could easily be accommodated by a SUSY theory with the phases somehow reduced 
by just the right amount. These ideas are clearly in conflict: one cannot have a theory 
which avoids fine-tunings by setting the SUSY parameters real, and at the same time 
expect dn near its current upper bound. The purpose of this paper is to describe a 
mechanism by which a moderate amount of SUSY CP violation can naturally appear 
in a theory in which the soft SUSY breaking terms have been taken real. 

The superpotential of the Minimal Supersymmetric Standard Model (MSSM) con- 
tains the Yukawa sector of the theory, Wy-, and a Higgs mixing term, 

Wmssm = Wy + /iH,Hd, (1) 

where H„ and Yld are Higgs doublet superfields. If the soft breaking terms come from 
the superpotential, as in (Q), then one can use H„ and Yid to rotate away the phase 
of /i. 

In order to avoid an additional hierarchy problem brought on by h/Mgut 
1 the MSSM is often extended by adding a singlet superfield N, whose scalar 
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component's vacuum expectation value (VEV) generates the Higgs mixing term (see 
and references therein). We refer to this model as the N+MSSM. One can use an 
R symmetry to forbid B and L violating terms in Wy-, and to allow only cubic terms 
involving N, so that the superpotential can be written as 

Wn+mssm = Wy + /iNH„Hd + aNl (2) 

Note that we can use the Higgs and singlet N superfields to rotate away the phases of 
h and a. This again assumes that the soft SUSY breaking Lagrangian can be written 
in the low energy supergravity (SUGRA) parametrization ||1CI|| : 



soft 



I |2 I |2 
\ipi\ 



+ B 



+ h.c. 



(3) 



where ipi are the scalar superpartners, A are the gauginos, and [ ]^ means take the 
scalar part. Here W^'^^ and W^^^ are the quadratic and cubic pieces of the superpo- 
tential, so that in the MSSM, W^^^ = Wy, and W^"^^ = /iH^H^; and in the N+MSSM, 
_ have defined the soft breaking parameters A = Attiq and B = Bra^ 

to include a mass scale mg. The parameters A, B, their mass scale mo, and the 
gaugino masses rhx, can all be complex. These parameters contribute to dn at the 
order of 10~^^(^/M^e cm, where is a combination of the phases of the parameters, 
and is a combination of superpartner masses, normalized to the weak scale. The 
only known ways to make such a large c?„ compatible with the experimental upper 
bound are to fine-tune the phase (f to order 10~^; have superpartner masses of order 



a few TeV; or somehow require all the phases to naturally be zero ||TT[. Both the first 
and second approach eliminate much of the attractiveness of SUSY [|1^. For exam- 
ple, having large superpartner masses virtually eliminates the possibility of radiative 



breaking of SU2 x f/i, which was one of the major successes of SUSY. Losing this is 
especially undesirable now that the top mass is large enough to make it work. So 
we will henceforth assume that A, B, tuq and fhx are all real [0. We do not have 



an explanation for how these conditions will be satisfied, but merely state that any 
complete SUSY model which has superpartners of order the weak scale must either 
satisfy these criteria, or provide an explanation for how their phases could naturally 
be of order 10~^ ill . 



Imposing these 'no fine-tuning criteria' means that the only source of CPV in 
either the MSSM or the N+MSSM is the CKM phase g]. CKM contributions to dn 
and de from renormalization group running |^ and from finite effects |TT| are below 
10~^° e cm, and are thus unobservably small. So if dn or dg were detected in the near 
future, could a SUSY theory with superpartners of order the weak scale explain them 
without resorting to fine-tunings? 

There have also been some interesting models of baryogenesis at the electroweak 



scale |1T5| , |T6|, |T^, which require CPV beyond the CKM phase [T^. Could one construct 
a model with sufficient CPV for electroweak baryogenesis, while satisfying the upper 
bounds on dn and de, without fine-tunings? 

With these questions in mind, we describe a mechanism by which a moderate 
amount of CPV can naturally arise in a non-minimal SUSY theory through loop 
corrections to the Higgs potential. The idea is that a phase which is unobservable 
at tree level can introduce an observable effective phase through loop effects. This 
effective phase will always be smaller than a tree level observable phase because of the 
usual factors of suppression associated with loops. Such a phase can make moderate 
contributions to dn and de, and may be useful in explaining the observed baryon 
asymmetry. 



In Section II, we present an illustrative model which provides an existence proof 
for this moderate CPV mechanism. In that section, we show that at tree level in the 
scalar potential, SUSY CPV is essentially unobservable. In Section III, we show how 
this CPV can appear in the observable sector in the one loop effective potential. The 
magnitude of this CPV will be suppressed by loop coefficients, so that dn and can 
naturally be near their current experimental bounds. 

2 The Model 

Let us construct a model which has complex couplings only to terms which contain 
singlet scalar fields which have zero VEVs. We also need these particles to have no tree 
level couplings to quarks or leptons. This means that the Higgs scalar potential will be 
CP conserving, as will all tree level vertices outside the neutral Higgs sector. At this 
order, there is no one loop contribution to d„ or dg. After one loop corrections to the 
scalar potential (V), a small phase can be induced into these vertices, which generates 
moderate dn and de- To do this in a model which is technically natural, one needs 
to add at least two such singlets {N', N") to the N+MSSM. In order that they have 
zero VEVs, we impose a discrete symmetry on their superfields: (N', N") — > — (N', N"). 
Then the most general cubic superpotential respecting this additional symmetry is: 

W = Wy + hNRJld + aN^ + c'NN'^ + c"NN"^ + 6NN'N". (4) 

One sees that the fields N' and N" have no direct couphngs to quarks, leptons, or 
gauge particles. We will see below that they can each have zero VEVs, and thus 
their couplings will not affect the tree level minimum of the scalar potential. They do 
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not affect CP violating observables studied to date (at tree level in V"), so we term 
this sector invisible. This is merely nomenclature. It should be possible to detect 
these particles, and perhaps even to see CP violating effects directly in processes in 
which they are produced, but they are certainly invisible when considering one loop 
processes involving only external quarks, leptons, and gauge bosons. 

Notice that we do not have enough freedom to rotate away all the phases in 
(^, and that after making the visible sector CP conserving, the reparametrization 
invariant h'^c'*c"* can be complex [l^. This is the phase which will produce CPV in 
the one loop scalar potential. Our first task is to be sure that this phase does not 
produce any CPV in the visible sector at tree level in V , else dn will again be too 
large. 

All supersymmetric contributions to dn come from the mass matrices of squarks 
and gauginos — if the mass matrices can all be made real, the SUSY contribution to dn 
disappears. If they are complex, the gaugino-squark-quark couplings become complex 
and contribute to dn through loop diagrams [|l|. Let us write the down squark mass 
matrix in a partially diagonalized basis: 

lidL^l + Ml {A* - /ine^^i tan (3)Md 

{A* - hne'^' tan PYMd f^dRl + Ml 

where M^, is the diagonal, real, Np x Np quark mass matrix (where Np is the number 

2 

of families), and ^iqL,R^ ~ ^^3/2 • Here n = \ {N)\ (so that hn takes the place of fi of 
the MSSM), and tan/5 is the ratio of Higgs VEVs. The angle 61 is one of the relative 
phases between the three VEVs, and is defined in (p^. If (as we have assumed) the 
soft breaking parameters are real, and the minimum of the scalar potential ^ is CP 
conserving, then this matrix is real. 



(5) 
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Next we can write the chargino mass matrix, M^+, 

(6) 



mw g2V2 
g2Vi hne^^^ 



in the basis of ||20| (with the argument of the Higgs VEVs rotated into ^i). Here rhw 
is the SU2 soft breaking gaugino mass, and g2 is the SU2 couphng constant. Again, 
if the minimum of is CP conserving, then this matrix is real. 

Since we have added three neutral fields to the MSSM (or two to the N+MSSM), 
the neutralino mass matrix, M^o, is 7 x 7. We extend the basis of PD[ with ipj^, ipN', 
and ipN"'- 
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where gi is the Ui coupling constant, and is the Ui gaugino mass. The angle ^3 
is also defined in ([T0|). The cross terms X and Y, which mix ipN' and ipN" with the 
visible sector, are proportional to the VEVs of A^' and A^", so that if these VEVs are 
zero, ipN' and ipN" decouple from the visible sector. The resulting 5x5 visible sector 
matrix is real, if sin6'i = sin 6*3 = 0. In that case, all of the mass matrices are real, 
and there is no new SUSY contribution to dn- 

To see if this is the case, we must consider the scalar potential V. We define two 
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Higgs doublets of the same hypercharge, and their VEVs, 



and the VEVs of the singlet fields 



' ' 

1^1 . 



(02) = (Hu) 



' ^ 



(8) 



(9) 



It turns out in this model that if n! = n" = 0, there are only three combinations of 
these VEV phases, 



(10) 



which appear in the tree level scalar potential. Elsewhere, a more general linear 
combination of 6i and ^3 (with integer coefficients) can appear. 
Let us write the scalar potential for our model. 



V = {h' \Nf + ml) |0i|' + (/i^ \Nf + ml) |02|' 

- '[AhN + 3ahN*^ + hc'*N'*^ + hc"*N"*^) {(j)\<p2) + H.c. 
+Ai(0l0i)2 + X^i^h)^ + A3(0l0i)(402) + X4{4h){4h) 
+ml \Nf + ml \N'f + ml \N"f + 9a^ \N\^ + |cf \N'\^ + |c'f |A^"|' + (11) 
+ (4 |cf + |6|') lATA^f + (4 |c'f + |6|') lA^A^'f 



+ 



AaN^ + AcNN'^ + Ac"NN"^ + AbNN'N'' 



+2{b*c' + be"*) \N\^ N'N"* + 3ac'N*''N'' + 3ac"N*''N"' + c'c"*N''N"*' + H.c 



I ]\T*2 -\Tt2 



where El 



8 



Ai = A2 = {gl + gl)/%, A3 = {gl - g\)l^. K = h^ - gl/2, (12) 

which is just the scalar potential for the N+MSSM |^ plus terms which involve N' 
and A^". The minimum of V can be written as 



{V) = {Vn+mssm) + K^m'^ + Ki^n'n" + K^^n"^ 

+Kion'* + K22nV + Ko^n"\ (13) 

where Kij depend upon all the other parameters. One can show that for any choice 
of the parameters in Vn+mssm, there exists a set of {b, c', c"} such that n' = n" = 
is a true minimum. We will assume that this condition is satisfied, so that (V) = 
{Vn+mssm)- 

Finally we must be sure there is no problem with spontaneous CPV. As we said, 
the potential depends only upon the three angles 6'i_3 (only two of which are inde- 
pendent). We can write 



{V) = ao — ai cos^^i — 0^2 cos^^2 — Q^s cos6'3, (14) 
where the are functions of the magnitudes of the three VEVs. Differentiating with 



respect to 61 and 6*3, we see that one solution to {V) = is sin6'j = 0. Romao p2 



showed that for this potential (i.e. {Vn+mssm)), this is the only stable minimum — 
that the spontaneous CP violating solution is actually a saddle point. Babu and 
Barr [^] make the interesting claim that this can be made into a minimum by large 
radiative corrections to the Higgs mass matrix, but these require very heavy squark 
masses (in which case hard CPV need not be suppressed by fine-tuning 0), small 
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charged Higgs mass, and tan (3 ~ These conditions make satisfying the CLEO 

bound on 6 ^ S7 nearly impossible It is also unlikely that a model satisfying 
these conditions could be consistent with such things as Grand Unification and solu- 
tions to the Dark Matter problem ||12|. Anyway, we can certainly choose parameters 
such that the minimum of V is CP conserving, and such that n' = n" = 0, so that 
all the SUSY mass matrices (H)-® are real at tree level. 

3 A Loop Induced Observable Phase 

It would seem that since the tree level potential is CP conserving, one could not have 
a one loop potential which is CP violating. The important point to remember is that 
even though the visible sector has no CPV, there are still CP violating couplings to 
A^' and A^". Consider, for example. Figure 1, which gives a purely finite contribution 
to a new term in V, 5A5(0|02)^ + H.c (this term is not present in the tree level 
potential, so it must be finite). The vertices are proportional to c'* and c"*, and the 
mixing between A^' and A^" contains pieces proportional to b. Thus 6X5 ~ b'^c'*c"*, 
which has a reparametrization invariant phase. For b, c', c" of order 1/2, 6X5 can be 
of order 10~^. 

Actually, the operator (to which Figure 1 contributes) is more accurately written 
as k {Nf (0l02)', which gives a contribution to {6V) of 2\k\ n^t>^f| cos(26'i + Arg/c), 
where Arg/c is just 6cp = Arg(6^ c'*c"*). We can write the general correction to {V) 
as 

integers 

{SV) = f^x,y,zCos{x9i + y93 + z9cp), (15) 

x,y,z 
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where the Kx,y,z are real coefficients, with the subscripts x, y, z taking on all integral 
values, though the Hx,y,z become negligible for large integers. Since our model is 
renormalizable p5| and {V) is CP conserving at tree level, all one loop terms in ([T5|) 



with z must be finite. Note that Figure 1 gives a finite contribution to (|1^) with 
{x,y,z) = (2,0,1). 

The perturbation in (|15D means that sin = is no longer a solution to {V{6i, 6^)') 
0. Since the Kx,y,z are small, the solution will lie close to this, so we can define 
6i = 0io + ^i) where sin^jg = 0- The minimization condition can then be written in 
terms of the hessian, and the perturbation: 



x,y,z 

^^V/^9^^9, d'V/d9l l\e,l V E^J:r' sin(x^io + l/^3o + ^^cp) 



(16) 

and solved for Ei. Note that sin(x6'io + y9^Q + z9cp) = ± sin(2;^^cp). To find the 
effective CP violating coefficient, recall that dn gets a contribution from the imaginary 
part of left-right squark mixing, which goes as sin^i ~ ^i. There will be finite 
contributions to from several terms in 6V, but they will be of the same order or 
smaller than that of 6X5 from Figure 1. From (|16D , one finds that Figure 1 gives 
El ~ I^AsI ^ sin 2/5 sin ^cp- The squark mixing, gluino mediated contribution to dn 



[O] due to this Ei can be written as 



ci„-10"^^ecm . {SX^l sm9cp, (17) 



V M ; ' ^' Ami 
where we have defined a SUSY mass scale 



^2 ^ _ ^18) 

mo rUg 
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There is a similar neutralino mediated contribution to which is suppressed by 
me/md, and a^j/as, but enhanced by the fact that sleptons tend to be hghter than 
squarks. If we take A = 1, colored superpartners ~ 300GeV, sleptons ~ ISOGeV, 
and all other superpartners ~ lOOGeV, one can have 



dn^ 10-^^ sin 9cpe cm, (19) 
4 ~ 10"^^ sin^cpecm. (20) 

These estimates depend upon the parameters and the mass scales in the theory, but 
the point is that the contributions entering at one loop are naturally much smaller 
than those from SUSY phases which contribute through tree level vertices. 



4 Concluding Remarks 

We have considered supersymmetric models which avoid excessively large contribu- 
tions to dn and d^ by requiring the 'no fine-tuning criteria' to be satisfied, i.e. that 
A, B, mo and rhx must be real fill. We showed that it is possible for moderate CP 



violating effects to be induced at one loop in models which have singlets with zero 
VEVs. We used an illustrative model with superpotential (^ and found that a stable 
minimum exists at {N') = {N") = 0. This means that using the tree level scalar 
potential, no CP violating effects would be detectable in conventional CP violating 
observables because all of the SUSY mass matrices are real. We demonstrated that 
this model introduces small CP violating phases into the one loop effective potential, 
so that one is left with a moderate contribution to dn and d^. One could easily have 
dn and d^ near their current experimental bounds in such a model, without the need 
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for fine-tuning or large superpartner mass scales. 



Note that SUSY contributions to dn from three gluon operators do not affect 
our conclusions. Assuming that the no fine tuning criteria are satisfied, such operators 
will also give a negligible contribution to dn- After one loop corrections to V in our 
model, there will be small contributions to dn from these operators, but they will 
probably be smaller in magnitude than the quark EDM contribution |Q. Thus ([T9|) 
and ( pO[ ) are reasonable estimates of the natural size of SUSY CPV possible in a 
model such as ours. 

We have discussed the issue of spontaneous CPV in Section II in the context of 
our model and concluded that we can easily choose the minimum of to be CP 



conserving. It is worth noting that Maekawa ||23| considered generating spontaneous 



CPV at one loop in the MSSM, though Pomarol ||2^ showed that such a model requires 
a CP odd Higgs which is too light. Pomarol also made the interesting point that a 
N+MSSM model (which does not rule out H„Hd, N, or terms by a symmetry) with 
a strictly CP conserving Lagrangian might violate CP spontaneously at tree level 
with a phase of order 10~^, and might be able to explain the e parameter as well as 
give dn near the current experimental bound [^. The trouble is that the fine-tuning 



needed by such a model of spontaneous CPV is actually much worse than that for 
hard CPV because the condition which must be satisfied is of the form 



cos^ 



X 
Y 



1 - (21) 



where 6 (or tt — ^) is the relevant spontaneous CPV phase, and X and Y are some 
combination of parameters and VEVs. We need 6 to be small to satisfy the bound 
on dn, which we can achieve only if 5 = (F — X)/Y is of order 6"^. For example, if 
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we need 9 ~ 10^^, then 5 must be fine-tuned to be of order 10"'', which is completely 
unacceptable. 

As we alluded to in the Introduction, having a moderate amount of CPV is neces- 



sary in models which generate the baryon asymmetry at the electroweak scale |jT8|. A 
recent interesting model of electroweak baryogenesis used CPV from the Higgs scalar 
mixing coefficient /if 2 EH) which can be defined as the coefficient of the (^(j)\4>2^ term 
in {V). It was pointed out [|T6| that /if 2 can be rotated out of the Higgs potential, 
but the resulting phase which appears in the gaugino mass matrices was then used by 



p!7[] . They found that with the small phase allowed by the limit imposed by there 



is probably sufficient CPV for the observed baryon asymmetry. Our results change 
these conclusions in two ways. At tree level, there is no phase in the gaugino mass 
matrices (after imposing the no fine-tuning criteria), and no way for Arg/if2 to cause 
CPV. Then, using the one loop effective potential in a model such as ours, there can 
be an effective phase £1 (= 61 — 9iq) introduced into the gaugino mass matrix of order 
El ~ lO^^Ocp- Although this is a large suppression, 9cp can be of order unity and 
so El should generate the same level of CPV as the phase used in [jl^, which was 
bounded by dn anyway |^ . 

From the standpoint of explaining the baryon asymmetry at the electroweak scale, 
or dn and de near their current experimental bounds, a loop induced observable phase 
provides an attractive alternative to the fine-tuning needed in the MSSM. If dn or df, 
were observed in the near future, and if superpartners were determined to be of order 
the weak scale, SUSY model builders would have to appeal to a mechanism such as 
ours, which naturally explains small effective SUSY phases. 
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FIGURE CAPTION 
Fig 1: One loop contribution to the operator 5X5{(j)\(j)2)'^. The 'X' indicates N'-N' 
mixing, which is required if SX^ is to contain a reparametrization invariant phase. 
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